Interaction of copernicium with gold: Assessment of applicability of simple density 

functional theories 
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Interactions of Cn (element 112) atoms with small Au clusters are studied using accurate ab 
initio scalar relativistic coupled cluster method for correlation treatment and two-component rel- 
ativistic density functional theory (RDFT) to account for spin-dependent relativistic effect. The 
results demonstrate the failure of RDFT with simple generalized-gradient and hybrid functionals in 
describing Cn-Au bonds in complex systems. 



Successful identification of superheavy element 112 
(Cn) by thermochromatography on gold surface [H, 0| 
had resulted in an avalanche of theoretical studies of the 
nature and properties of the interactions of the Cn atom 
with gold (see and references therein). It is gen- 

erally believed (see, e.g., [1]) that the relativistic density 
functional theory (RDFT) with simple generalized gradi- 
ent approximations (GGA) for the exchange-correlation 
functional (XCF) ^ ij} provides a satisfactory tool for 
theoretical modeling of Cn/Au surface complexes. Since 
the quantitative experimental information on Cn - Au in- 
teractions is restricted to a single measured value, the ad- 
sorption temperature of Cn on gold surface, thus being 
obviously insufficient to estimate the reliability of theo- 
retical approaches, the main argumentation for the valid- 
ity of RDFT/GGA in applications to their description re- 
lies on the data concerning diatomic systems (see e.g. [1]) 
and on the experience of calculations of other heavy el- 
ement compounds. It should be noticed, however, that 
the bulk of such experience is nearly useless because of 
unique features of the Cn atomic structure in the valence 
region. The presence of relatively compact filled d-shell 
with rather high orbital energies (resembling that in Au 
and Hg) enables one to expect unusually strong d^'^ — d^^ 
contributions to the Cn-Au bonding, in a sense similar to 
the aurophilic attraction [gj. The efficiency of the RDFT 
with semilocal functionals in such cases is at least ques- 
tionable as demonstrates the example of the Au„ cohe- 
sion energy which is nicely reproduced by RDFT/GGA 
in the Au diatomic but becomes progressively underesti- 
mated for larger clusters and bulk [13, [HI . The depen- 
dence of the accuracy of RDFT/GGA based results for 
gold clusters on the coordination numbers of Au atoms 
has been shown in Ref. [lo| . 

The mentioned facts inspire serious doubts concerning 
the adequacy of the RDFT/GGA treatment of the sys- 
tems under discussions. These doubts are strengthened 
by the detailed pilot comparison of the results of relativis- 
tic DFT and many-body perturbation theory (MBPT) 
calculations on moderate-size Hg-Au„ and Cn-Au„ com- 
plexes . Unfortunately, finite-order MBPT binding en- 
ergies for similar (Au- and Hg-containing) systems are 



strongly affected by convergency problems and the re- 
summation scheme chosen in this work was not optimal 
121. The errors due to the use of rather restricted one- 



electron basis sets with insufficient flexibility of high an- 
gular momentum components can also be non-negligible. 
A more accurate correlation treatment is, therefore, re- 
quired to pronounce the decisive verdict. 

In the present Letter we report the calculations on 
Cn-Au„ complexes, n = 1—4, by a hybrid scheme com- 
bining ab initio scalar relativistic correlation calcula- 
tions with the estimation of spin-dependent relativistic 
effects (effective spin-orbit couplings) through geometry- 
dependent corrections to interaction energies (A50) ob- 
tained at the RDFT level Q . The necessary accuracy of 
the correlation treatment is ensured by the use of an ac- 
curate coupled cluster technique and the extrapolation to 
the complete basis set limit. We have to stress that our 
scheme does not imply the smallness of spin-dependent 
interactions or the neglect of their interference with elec- 
tron correlations [l^ . The geometries of polyatomic sys- 
tems {n = 2—4) were restricted to those resembling pos- 
sible structures of adsorption complexes on the gold (111) 
surface, namely, the C2V configurations of Cn-Au2 and 
Cn-Au4, as well as the C^v configuration of Cn~Au3 sim- 
ulating the "bridge" and "hollow" positions of the Cn 
adatom, respectively. The Au-Au distances were fixed 
at their bulk value (2.884 A); only the distances r be- 
tween the Cn atom and the center of the Au subsystem 
were optimized. 

The computational scheme employed in the present 
study resembles that used in Ref. [ij]. The "small" 
atomic cores (60 electrons for Au and 92 electron for Cn) 
were replaced by accurate shape-consistent semilocal rel- 
ativistic pseudopotentials (RPPs) [15,, 161 optimized for 
the description of valence shells. Scalar relativistic calcu- 
lations were performed with the spin-averaged version of 
these RPPs by the conventional coupled cluster method 
with fully optimized single and double and perturbative 
triple amplitudes, CCSD(T), implemented in the CFOUR 
program package [l3|- For the systems with odd num- 
bers of gold atoms unrestricted Hartree-Fock references 
were used. The innermost explicitly treated atomic shells 
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{6s6p for Cn and 5s5p for Au) were not correlated. 

We used the correlation-consistent basis sets ([l3, El) 
constructed according to the prescriptions from Ref. [19j . 
Taking into account the experience of calculations on Au- 
Au and Au-Hg [l^, , the bases were augmented by dif- 
fuse functions (the "aug-cc" variant). The complete basis 
set limit (CBS) estimates were derived from the energies 
obtained with triple and quadruple zeta quality bases us- 
ing the conventional two-point formula 20] . As has been 
pointed out in Ref. [14j] , the extrapolation procedure effi- 
ciently suppresses basis set superposition errors, so that 
we had no need to invoke any counterpoise procedure. 

The contributions from the spin-dependent relativis- 
tic effects to the interaction energies were evaluated by 
comparing the results of one- and two-component (Ic and 
2c) non-collinear RDFT calculations performed with the 
code The data presented in this letter were ob- 

tained with simple generalized gradient approximation 

0] usually denoted by the acronym B88P86 and be- 
lieved to be well suited for superheavy element electronic 
structure calculations, and with PBEO hybrid functional 
[2^ chosen because of its (partially) non-empirical na- 
ture. 

The main results of the calculations are listed in Ta- 
ble m The equilibrium separations and binding energies 
derived from Ic CCSD(T)/CBS-|-Aso potential curves 
for the CnAu diatomic molecule are in very reason- 
able agreement with those obtained by RDFT with both 
XCFs employed. When the cluster size increases, the 2c 
RDFT binding energies deviate progressively from their 
Ic CCSD(T)-|-A5o counterparts; for the largest system 
under study the underestimation of the binding energy 
reaches ca. 40-45 %. It might be interesting to notice 
that the results of a similar numerical experiment with 
Cn replaced by the neighbor element 113 (El 13) with 
markedly lower energies of the filled d shell do not ex- 
hibit a similar trend. 

The extrapolated MBPT-I-A50 binding energies for 
Cn-Au4 from Ref. Q are visibly too large; neverthe- 
less, the main conclusion of the cited work concerning 
the failure of RDFT with simple GGA and GGA-based 
hybrid energy functionals in describing more or less com- 
plex Cn-Au systems is confirmed by the present accurate 
calculations. A reasonable agreement between the most 
recent RDFT/GGA Cn/Au adsorption energy estimate 
(0.46 eV, [3]) and the presently assumed experiment- 
based value (0.54^3 eV, [2|) seems not to be sufficient to 
prove the opposite opinion and, probably, results from a 
fortuitous cancellation of errors in ^Sj (for instance, those 
of approximate nature of the B88P86 functional and of 
use of purely local functional to get the electron and mag- 
netization density distributions). 

The work was supported by the Russian Foundation for 
Basic Research (Grant No. ll-03-12155-ofi-m-2011). 



TABLE I: Equilibrium r values (A) and Cn-Aun binding en- 
ergies (eV). 
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*^ Ref. ;5]; fourth-order MBPT was followed by the 
extrapolation using constrained [3/2] Pade approximant. 
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FIG. 1: The dependencies of CCSD(T)-HAso and RDFT 
binding energies on gold cluster size. The data on E113-Au„ 
systems are taken from Ref. p^ . 
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